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Abstract 

Starting from a small number of well- motivated axioms, we derive a 
unique definition of sums with a noninteger number of addends. These 
"fractional sums" have properties that generalize well-known classical sum 
identities in a natural way. We illustrate how fractional sums can be used 
to derive infinite sum and special functions identities; the corresponding 
proofs turn out to be particularly simple and intuitive. 

"God made the integers; all else is the work of man." 

Leopold Kronecker 

1 Introduction. 

Mathematics is the art of abstraction and generalization. Historically, "num- 
bers" were first natural numbers; then rational, negative, real, and complex 
numbers were introduced (in some order). Similarly, the concept of taking 
derivatives has been generalized from first, second, and higher order derivatives 
to "fractional calculus" of noninteger orders (see for instance [9]), and there is 
also some work on fractional iteration. 

However, when we add some number of terms, this number (of terms) is still 
generally considered a natural number: we can add two, seven, or possibly zero 
numbers, but what is the sum of the first —7 natural numbers, or the first tt 
terms of the harmonic series? 

In this note, we show that there is a very natural way of extending sum- 
mations to the case when the "number of terms" is real or even complex. One 
would think that this method should have been discovered at least two hundred 
years ago — and that is what we initially suspected as well. To our surprise, 
this method does not seem to have been investigated in the literature, or to be 
known by the experts, apart from sporadic remarks even in Euler's work [5] (see 
equation (fT2l) below) . Of course, one of the standard methods to introduce the 
r function is an example of a summatioja with a complex number of terms; we 
discuss this in Section FOl equation ([5]) 



^Note by the second author. Many of tlie original ideas in this text are due to Marlcus 
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Since this note is meant to be an introduction to an unusual way of adding, 
we skip some of the proofs and refer the reader instead to the more formal note 
0- Some of our results were initially announced in [7]. 

1.1 The Axioms. 

We start by giving natural conditions for summations with an arbitrary com- 
plex number of terms; here x, y, z, and s are complex numbers and / and g 
are complex-valued functions defined on C or subsets thereof, subject to some 
conditions that we specify later: 

(51) Continued Summation 

ly—x u—y^l v—x 

(52) Translation Invariance 

E /M = E/(^ + ^)- 

v—x-\-s U—X 

(53) Linearity for arbitrary constants A, ^ G C, 

y y y 

E (A/(^) + ^^9{^)) - A E /(^) + E 

y—x v—x v—x 

(54) Consistency with Classical Definition 

1 

E/M = /(!)■ 

v=l 

(55) Monomials for every d € N, the mapping 

z^E-' 

v=\ 

is holomorphic in C. 

Miiller, who invented them while he was a high school student in the remote German province 
town of Morsbrunn, Bavaria. He was lacking the skills to carry out a formal mathematical 
proof, but he kept producing the most obscure mathematical identities on classical sums and 
fractional sums. I met him at the science contest "Jugend forscht" for high school students, 
and from then on helped him to turn his ideas into actual mathematical theorems and proofs, 
and to find out which of his formulas and identities were correct (most of them were). The 
main results have been published in The Ramanujan Journal [8]. 
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(S6) Right Shift Continuity if lim„^_|_oo /(^ + n) = pointwise for every 
z e C, then 

y 

lim 5]/(i. + n)=0; (1) 

v—x 

more generally, if there is a sequence of polynomials (pri)neN of fixed degree 
such that, as n — > +cxd, |/(z + n) — p„(z + n)\ — > for all z G C, we 
require that 



0. 



(2) 



The first four axioms (S1)~(S4) are so obvious that it is hard to imagine 
any summation theory that violates these. They easily imply YTv=\ /(^) — 
/(I) + /(2) + • • • + /(n) for every n € N, so we are being consistent with the 
classical definition of summation. 

Axiom (S5) is motivated by the well-known formulas 



n(n + 1) 



n(n+ l)(2n+ 1) 
6 



E-^ 



n(n + 1) 



and similarly for higher powers; we shall show below that our axioms imply that 
all those formulas remain valid for arbitrary n G C. 

Finally, axiom (S6) is a natural condition also. The first case, in ([T]), ex- 
presses the view that if / tends to zero, then the summation "on the bounded 
domain" should do the same. In ([2]), the same holds, except an approxi- 

mating polynomial is added; compare the discussion after Proposition ll.il 

It will turn out that for a large class of functions /, there is a unique way 
to define a sum f with z G C that respects all these axioms. In the next 
section, we will derive this definition and denote such sums by J ■ We call 
them "fractional sums." 



1.2 Prom the Axioms to a Unique Definition. 

To see how these conditions determine a summation method uniquely, we start 
by summing up polynomials. The simplest such case is the sum X]^=i with 
c G C constant. If axiom (SI) is respected, then 

1/2 1 1 

E^+ E c = E^- 

v=l iy=3/2 1^=1 

Applying axioms (S2) on the left and (S4) on the right-hand side, one gets 

1/2 1/2 

Ec+E^ = c- 

i,= l u=l 
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It follows that — "^Z^- This simple calculation can be extended to cover 

every sum of polynomials with a rational number of terms. 

Proposition 1.1. For any polynomial p : C — )• C, let P : C ^ C be the unique 
polynomial with P(0) = and P{z) — P{z — 1) = p{z) for all z € C. Then: 

• The possible definition 

V 

Y,p{v):^P{y)-P{x-l) (3) 

u—x 

satisfies all axioms (SI) to {S6) for the case that f is a polynomial. 

• Conversely, every summation theory that satisfies axioms (SI), (S2), (S3), 
and (S4) also satisfies (0) for every polynomial p and all x,y € C with 
rational difference y — x d Q. 

• Every summation theory that satisfies (SI), (S2), (S3), (S4), and (S5) 
also satisfies for every polynomial p and all x,y € <C. 

Proof. To prove the first statement, suppose we use ([3]) as a definition. It is 
trivial to check that this definition satisfies (SI), (S3), (S4), and (S5). To see 
that it also satisfies (S2), consider a polynomial p and the unique corresponding 
polynomial P with P(a:;) — P(a;—f) = p{x) and P(0) = 0. Define ^(x) ■.= p{x + s) 
and P{x) := P{x + s) - P(s). Then P(0) = 0, and P{x) - P{x - I) = P(a; + 
s) — P{x + s — I) = p(x + .s) = p{x). Hence 

V y y+s 

^p(z. + s) = ^p(i.) = P(y)-P(x-I) = P(2/ + s)-P(x + s-I)= J2 P^")- 

1/—X v—x v—x-\-s 

To see that ([3]) also satisfies (S6), let V„ be the linear space of complex polyno- 
mials of degree less than or equal to ti e N. The definition |jp|| := Y^i=Q 
p & Va introduces a norm on Va-. If we define a linear operator : ^ C via 
Y^xP ■~ jyt=xP(^)^ then this operator is bounded since dimVo- — a + 1 < oo. 
Thus, if (g„)„gN C Vo- is a sequence of polynomials with lim„_^oo IknII = 0, we 

have Ihxin^oo 'Y^^.j. Qn{i^) — 0. Axiom (S6) then follows from considering the 
sequence of polynomials (gn)neN with g„(a;) := p{x + n) — p„(a; + n) and noting 
that pointwise convergence to zero implies convergence to zero in the norm || • || 
of (7„, and thus of 

To prove the second statement, we extend the idea that we used above to 

^ 1/2 

show that X)i/=i'' = c/2. Using (SI), we write for an integer r > 1 

r rjs 1r I s r 

E-' = E-'+ E + E 
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where the left-hand side has a classical interpretation, using (SI), (S2), and 
(S4). Rewriting the right-hand side according to (S2) and using (S3), we get 

r s-1 r/s , l, \ d ^/^ s-1 r/s 

v=l k=Ou=l ^ ^ iy=l fc=Oiy=l 

where the qd-i,k{v) ~ {y + kr/sY ~ i^"^ are polynomials of degree d — 1 (and 
all 9-1, fc = 0). Now we argue by induction. If rf = 0, the previous equation 
clearly determines X]I/=i 1 and by linearity also the corresponding sum over 
arbitrary constants c e C. Once the value of the sum of any polynomial of 
degree c? — 1 is determined, the equality also determines the value of J^Jji v'^^ 
and by linearity, the sum of every polynomial p of degree d. Using (S2) again, 
we see that the axioms (SI) to (S4) uniquely determine the sum X]1^=k-P(^) — 
jy^=i^^ p{'^ + X — 1) ii y — X e Q. As we have seen, equation ([3]) is a possible 
definition satisfying those axioms; hence it is the only possible definition for 
y-x€Q. 

Finally, it is clear how the restriction y — x^Q can be lifted by additionally 
assuming (S5) (equation (j3|) is already satisfied for x — y e R by requiring just 
continuity in (S5); holomorphy is required for x — y G C). □ 

Consider now an arbitrary function /: C — > C. If we are interested in 
Y^^=x /(^) ^'^^ complex x,y G C, we can write 

V y+n x+n—l v+n 

E/H+ E /H= E E /W' 

U—X u—y-]-! U—X u—x+n 

where n € N is an arbitrary natural number. Hence, 

y x+n-1 y+n y+n 

E/('^) = E /M- E /M+ E 

U—X U—X I/— y+1 f— x+n 

n y 
= ^(/(^. + x-l)-/(i. + y))+^/(i/-Hn). (4) 

U—l U—X 

What have we achieved by this elementary rearrangement? In the last line, the 
first sum on the right-hand side involves an integer number of terms, so this can 
be evaluated classically. All the problems sit in the last sum on the right-hand 
side. The payoff is that we have translated the domain of summation by n to 

the right. Since Q holds for every integer n, we can use (S6) to evaluate the 

limit as n — > oo: if /(n + z) — > as n — > oo for all z, then ( S6 ) implies that the 
limit as n — > oo of the last sum should vanish. We get 

V oo 

E/M = E(/(^+^-i)-/(^ + y))- 

U=X U—1 
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This is of course a special condition to impose on /, but the same idea can 
be generahzed. For example, if /(z/) = Inz/, then for v G [x, y\ C K"*", the values 
/{v + n) are approximated well by the constant function /(n), with an error 
that tends to as n oo: we say that / = In is "approximately constant." 
Using (S3), 

y y y 

hi{h' + n) — In n + ( ln(:^ + n) — In n) 



for every n e N. But by (S6), the last sum vanishes as n oo, while the 
first sum on the right-hand side has a constant summand and is evaluated using 
Proposition ll.il Taking the limit n — oo in (j4]), it follows by necessity that 

y / n \ 

Inv ~ lim j (ln(i/ + x — 1) — ln(i/ + y)^ + {y — x + l)lnn \ . 

v=x \z^— 1 / 

Before generalizing our definition further, we take courage by observing that 
this interpolates the factorial function in the classical way: we define 

ri/M :=exp(f]ln/M] 

U—X \u — x J 

and thus get 

z\ = TT V — lim exp ( In ( ) + z In n | 

11 „_>oo ^ I ^ \V + Z ] 



v=\ \v=\ 



1™ u^-n^Hr(^+i)' 



using a well-known product representation of the F function [1] 6.1.2]. 

It is now straightforward to use the heuristic calculation in Q together with 

Proposition 11.11 and axiom ( S6 ) to derive a general definition: all we need is 
that the value of /(n-|-z) can be approximated by some sequence of polynomials 
p„(ri -I- z) of fixed degree for rt — > oo. 

Some care is needed with the domains of definition: the example of the 
logarithm shows that it is inconvenient to restrict to functions which are defined 
on all of C. All we need is a domain of definition U with the property that z G U 
implies z -I- 1 S U. This leads to the following (using the convention that the 
zero polynomial is the unique polynomial of degree — oo). 

Definition 1.2 (Fractional Summable Functions). 

Let U C C and cr e N U {— oo}. A function f : U ^ C will be called fractional 
summable of degree <j if the following conditions are satisfied: 

• X + 1 eU ior SiW X eU] 
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• there exists a sequence of polynomials {pn)neN of fixed degree a such that 
for all X G U 

\f{n + x)—pn{n + x)\ — > as n — +c» ; 

• for every x,y + 1 G U, the limit 

n+y n \ 

^ p„(i.) + ^(/(i. + a:-l)-/(j. + y)) 

u=l J 

exists, where is defined as in ([3|l. 

In this case, we will use the notation 

y y 
g/(zy) or briefly ^/ 

L^ — X X 

for this limit. Moreover, we can define fractional products by 
^fii^) :=exp(^ln/H] , 

whenever In / is fractional summable. 

Note that this definition does not depend on the choice of the approximat- 
ing polynomials {pn)neN- if iPn)neN is another choice of approximating poly- 
nomials, then lim„^oo ipn{n + x) ~ Pn{n + x)) = for all x £ U, and hence 
for all a; G C since the set of polynomials of degree at most tr is a finite- 
dimensional linear space. As shown in Proposition 11.11 sums of polynomials 

satisfy axiom (S6). Substituting for / and Pn ~Pn for pn in (S6) proves that 

lim„^oo [E'Zl+xPni^) - E]^ll+xPni'')) = 0- 

Moreover, this definition is the unique definition that satisfies axioms (SI) 
to (S6): 

Theorem 1.3. Definition ll.Si satisfies all the axioms (SI) to {S6) (for suitable 
domains of definition), and it is the unique definition with this property (for the 
class of functions that we are considering). 

Proof. We have already proved uniqueness above, by deriving Definition 11.21 

from the axioms (SI) to (S6). It remains to prove that this definition indeed 
satisfies all the axioms. Clearly, (S3) and (S5) are automatically satisfied. Sub- 
stituting the definition into (SI), (S2), and (S4), these axioms can be confirmed 

by a few lines of direct calculation. To prove (S6), we use the definition and 
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the other axioms (in particular continued summation (SI)) and calculate 

Cy ^ \ 

^ f{iy + 72) - Y^Pnif + n)] 
y—x U—X / 

(y y n \ 

E /(i' + n) - E f{u) + E {f{^ + x-l)~ f{v + y)) 
y—x u—x v—l / 

(y+n y n+a;-l y+n \ 

E /(-) - E + E - E = 0- 
u—x+n u—x v—x v—y+\ J 

This proves that Definition 1 1 . 2 1 satisfies all the axioms. □ 

2 Properties of Fractional Sums. 

Now that we have a definition of sums with noninteger numbers of terms, it is 
interesting to find out how many of the properties of classical finite sums remain 
valid in this more general setting, and what new properties arise that are not 
visible in the classical case. 

2.1 Generalized Classical Properties. 

One of the most basic identities for finite sums is the geometric series. For 
simplicity, let < g < 1. Then the function v ^ ys, approximately zero (we 
have lim„_>oo 5^+" = for every 2 e C), and the definition reads 

a; 00 00 x+\ 

E '^'^ = E ('^"'^ - = (1 - 1'^') E ^"-^ = -T^- (6) 

i/=0 u=l i/=l ^ 

Thus, the formula for the geometric series remains valid for every x G C. 

A similar calculation shows that the binomial series remains valid in the 
fractional case: for every c G C \ { — 1, —2, —3, . . .} and x £ C with |a;| < 1, we 
have 

There are generalizations of ([5]) to the case q > 1 and of (O to the case 
|a;| > 1: these involve a "left sum" as introduced in Section |3l 

An example of a summation identity with more complicated structure is 
given by the series multiplication formula 

(e/h) • (e^m) 

x / u — 1 u — 1 \ 

= E /MffM + zMEffW+sME/w (8) 

v=l \ fc=l fc=l / 
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for every x G C, given that all the three fractional sums exist (see [SI Lemma 7]; 
it generalizes the formula (oi + a2)(6i + 62) — ai&i + 02^2 + a2&i + 6201, and 
similarly for all positive integers x. 



2.2 New Properties and Special Functions. 



As shown in Section [1.2| our definition interpolates the factorial by the F func- 
tion, 



z! EE g = r(z + 1). (9) 
i/=i 



An amusing consequence is 



-1/2 

^ (i^^ + 1) = tanhTr; (10) 



this is because 



-n\ ^ ,N -TTr -n\ ■^ r(i/2 + i)r(i/2 - z) 
+ - H(- + ^)H(--^)= r(i + .)r(i-.) 



r(l/2 + i)r(l/2-i) sin(7ri) sin m 



ir(i)r(l — i) i sin(7r(l/2 + i)) icosiir 

sinh TT 

tanhTT 



coshTT 



using r(z)r(l — z) = 7r/sin(7rz) = in/ smhlniz). 

Many basic fractional sums are related to special functions. As a first ex- 
ample, consider the harmonic series. Since 1-> i^^^ is approximately zero, the 
definition reads 



and in particular 



E- = -2 + + =-21n2, (12) 



V V 2 3 4 



which was noticed already by Euler [5l pp. 88-119, §19]. For general x, the 
harmonic series can be expressed in terms of the so-called digamma function [1] 
6.3.1] i)(x + \) = ^lnr(2: + l) and the Euler-Mascheroni constant 7 = 0.577. . .: 
one obtains [IJ 6.3.16] 

X °° /I 1 \ 



1^ '■ — ' \ V V + X 
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Note that the reflection formula T, 6.3.7] for the digamma function becomes 

^ - = 7rcot(7ra;). (14) 



As a further generalization, it is convenient to consider the Hurwitz C, func- 
tion, traditionally defined by the series 



oo _^ 

C{s,x):=y — 



5R(s) > 1. 



By analytic continuation, C(s,x) can be defined for every s S C, except for a 
pole at s = 1. For a: = 1, the Hurwitz Q function equals the well-known Riemann 
^ function: 

C(s,i) = C(s). 

It turns out that the Hurwitz C, function can be understood as a fractional power 
sum. It can be shown 8, Corollary 14] that for every a G C \ {—1} and for all 
a; eC\ {-1,-2,-3,..'.}, 



^Q{-a)-C{-a,x + l). 



(15) 



A useful special case is 



^v'^ = {2- 2-'^) C(- 



(16) 



Note that such equations give in many cases intuitive ways to compute properties 
and special values of special functions. Everybody knows the formula X^^^j^ v = 
x{x + l)/2, so -E^if V = -1/8, and thus by ([T6l) 

-1/2 



i=£.^ = (2-2-)C(-l) = |c(-l) 



It follows that C(— 1) = —1/12. Similarly, we have 

-1/2 

=0 dz 



^(2-2-)C(-.) 



-1/2 

^ Ini^ = In ^ 1/ = InT - 



(in the second equality, we interchanged differentiation and fractional summa- 
tion; it is not hard to check that this is indeed allowed). Since r(l/2) = ^Jti, 
this easily implies that C'(0) = -(1/2) ln(27r). 
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Similarly, differentiating (jTS]) b times with respect to a and arguing as before 
(compare also [H Sec. 6]), we obtain 

X 

g z.'^(lnt.)^ = (-1)'' [c^'\-a) - C^'H-a,x + 1)) . (17) 
There are some classically unexpected special values like 

E-i^^- = -^-^c'(-i). (18) 

1^—1 

2.3 Mirror Series and Left Summation. 

There is an identity for classical sums which is almost never mentioned, because 
it seems so trivial. Consider the sum 

/(-10) + /(-9) + /(-8) + /(-7). 

Obviously, there are two formally correct possibilities to write this sum, 

-7 10 

either ^ f{iy) or ^fi-i")- 

Classically, it is clear that Y^v=a — J2^=-b /(~^)- Does this carry over to 
the fractional case? There is a fundamental problem: our definition of fractional 
sums involves lim„_>+oo f(v + n), i.e., / is evaluated near +oo, and when f(v) 
is replaced by f{—v) then / would be evaluated near —oo where the values may 
be unrelated. This will be discussed in the next section. 



3 An Alternative Axiom and Left Summation. 

Looking back at the axioms given in Section II. 1) there is one axiom that could 

possibly be modified: in (S6), limits as n — > +oo are considered, but one could 
equally well look at limits as n — >■ — oo. This way, one obtains an axiom of "left 
shift continuity" : 

(S6) Left Shift Continuity if lim„^oo f{z~n) — pointwise for z e C, then 



lim y f{h' — n) = 0; 



(19) 



more generally, if there is a sequence of polynomials {pn)nm of fixed degree 
such that \f{z — n) — p„(z — n)| — > for z g C as n — >■ oo, we require 
that 

y y 
^/(i^-n)-^p„(zy-n) — >Q. 
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Repeating the calculations of Section 11.21 one gets an alternative definitioiH 
which we do not state here formally: it is exactly the same as Definition II. 2i 
except that in every limit, n — >■ oo is replaced by n — — oo. 

It can be shown that this definition is the unique one that satisfies axioms 

(SI), (S2), (S3), (S4), (S5), and (S6). Note that in general, the existence of g 
and ^ are independent, and if both left and right fractional sums exist, they 
may have different values. For example, for every z with 5R(2;) > 0, we have 

-1/2 

£ v'-^{-\r^^ (2-2-^)C(-z), 

in contrast to equation (jl6p . 

What we do have is the obvious relation 

E/M= E /(-^)- (20) 



4 Classical Infinite Sums, Products, and Limits. 

Fractional sums are not simply a new world with results that have no meaning 
in the classical context; they allow us to derive identities that can be stated 
entirely in classical terms. Some of these formulas are known and some seem 
to be new. Of course, all these identities can in principle be computed without 
fractional sums. But proving them with the help of fractional sums is rather 
intuitive and simple, since most of the steps use fractional generalizations of 
basic, very well-known classical summation properties. 



4.1 Some Infinite Products. 

As a first example, we show how to compute a closed-form expression for the 
infinite product 

P(x):= lim n 1 + T (21) 

k=l ^ ^ 

for X > —1/2. It was first considered by Borwein and Dykshoorn in 1993 (see 
[3]). By taking logarithms, one gets 



InF(x) ^ - E fsfcln (l + - 2 f fc - 1) In (l 



k=l 



2k J V 2/ V 2(fc-l/2) 



2x 



^Note that no other complex directed Hmit to infinity (like n — > ioo) can determine a 
definition uniquely: only adding or subtracting n G N to the upper summation boundary 
consists of adding or subtracting n terms to the series, which can be done classically. 
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Consider the function H> 2i/ln (1 + xjv) that tends to 2a: as — ?> oo. Accord- 
ing to Definition 11.21 we have 



-1/2 

^ 2^ In [l 



hm 



1.2. + |:2H.(l + i)-2(.-i),„(l + ^) 



Thus, we get 
lnP(a;) 



2 ^ j/ln 



-1/2 -1/2 

E2.1n(l + ^)=- 

-1/2 -1/2 

2 ^ J/ ln(z^ + + 2 ^ v\nv 



V 



-1/2+2; 



In 2 



2 £ (^-.^)ln^- 3C'(-1) 



1^—1+2 
-1/2+2; 
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-1/2+2; 



In 2 



2 ^ iylni/-|-2a; ^ In - — - 3C'(-1) 

I/— l+X 



1^=1+2; 

"2:-2 ( C' f-l,x + 



2x In \x 



1 



ln(a;! 



In 2 
T2" 



3C'(-1), 



where we have used equation p^ . index shifting, equation ([T71) . equation JS]), 
and continued summation. By exponentiating, we finaUy get 



22; 



P 



ix) = 2"T2 I LI^ILII I g-2:-2C'(-l,2;+i)+2C'(-l,2;+l)-3C'(-l) ("22) 

I r(a; + 1) / 



Using Mathematica's built-in numerical procedures, this infinite product 
identity can be checked numerically. Figure [1] shows a comparison of both 
sides of this equation. 

By application of this method, a large class of infinite products can be ex- 
plicitly computed, which seems to include the class of products considered in [2]. 
Here is an example of a new identity: using the same steps as in the calculation 
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Figure 1: Numerical check ol (|22p . The right-hand side corresponds to the 
lowest curve, while the other three curves (from top to bottom) are plots of 
approximations to P{x) (finite products as in (HIJ) for n = 1, n = 10, and 
n — 50 respectively. 

above, one easily proves that for x > — 1, 



-1/2 



^2kln'{2k + x) = 41n2(c'(-l,^]-C'fl,^ + l 



k=l 



X 

-^'2+1 



r 

_2<;"f_l,i±IUx<"(o,£±i 



hi 2 



<"(o,| + i) 



Resolving the definition and exponentiating, we get the following classical 
limit identity: 



lim 

n— »oo 



(2n)"5-=^-(«+3)M2«) Y\{k + a;)(-i)'*=Mfe+:r 



fe=i 



-2x In 2 



2-i ln2+4C'(-l,^)-4C'(-l,f + l) 



r(f + 1) 



g2C"(-i,f+i)-2C"(-i,^)+x(c"(o,5^)-C"(o,f+i))^ (23) 



Again, we have used Mathematica for a quick numerical check that is shown in 
Figure [H 
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Figure 2: Numerical check of (|23p . The right-hand side corresponds to the 
uppermost curve, while the other three curves (from bottom to top) are plots 
of the left-hand side for n = 10, n = 100, and n = 1000 respectively. 



4.2 The Multiple F Function and Series Multiplication. 

In this section, we consider the multiple gamma function r„, a generalization of 
the classical gamma function F, defined for n e N and z e C by the recurrence 
formula (compare 2 ) 



-+'^' + ^) - T;:(;r' 

Fo(z) = z-\ 
F„(l) = 1. 



(24) 



These equations do not determine the functions F„ uniquely, so one needs the 
additional Bohr-MoUerup-like condition that F„(x) is positive and n times dif- 
ferentiable on x e IR+, and that (^l)"^"'^^^ liiF„(a;) is increasing (see [4]). For 
n = 1, this definition reproduces the classical gamma function: Ti[z) = T(z). 
The (reciprocal of the) special case n = 2 is known as the Barnes G function 



By ([M]), it satisfies 



G{z) = 1/F2(z). 



G{n) =F(1)F(2)F(3)---F(n-1). 



More generally. 



l/F„+i(z) = F„(1)F„(2)F„(3) • • • F„(z - 1), 



so F„_|_i is the reciprocal of the product of F„, which means that F„(z) is 
something like an n-fold product of the first z — I natural numbers: 



(r„(^)) 



n n n n 

2^1— 1 1/2 — 1 ^^3 — 1 l^Ti — 1 
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While this equation only makes sense for z G N \ {0}, one can easily show that 
the definition of r„(z) for z G C is compatible with our definition for fractional 
sums and products, i.e., that for every 2; G C (except for poles) and n € N\ {0}, 
we have 

2 — 1 1^1 — 1 1^2 — 1 — 1 

(~ir+iinr„(z) = g E E ••■ E 1"^- 

We will now show some properties of the multiple gamma function r„, specif- 
ically for the example n — 2, simply by using basic fractional sum identities, 
without using any special function properties of r„. By the multiplication for- 
mula ([5]), we have 

z-l / jy-1 

lnG{z) = -lnr2(z) = ^ I 1-^lnfc 

u=l \ fc=l 

CZ~1 \ /z-l \ 2-1 2-1 / 1^-1 

El El-- -E1-1---E 1--E1 
jy = l / = l / I/=l u = l \ k = l 

2-1 

= (z- l)lnr(z) -gi/lnt/ 
i/=i 

= (z-i)inr(z)-HC'(-i)-C'(~i,^)- 

The last equality follows from equation (|17p . Thus, we have found an explicit 
formula for G{z) in terms of derivatives of the Hurwitz ^ function. 
Equation (fTSl) gives the special value 

These are of course very well-known results, but the calculations are strikingly 
simple. Moreover, this example shows that there is a wide variety of interesting 
"special functions" that do not have to be defined separately, but can be treated 

in a unified manner by our theory of fractional sums. New generalizations 

2-2 

comparable to G{z) = ^ nl include 
-1/2 

^(2n)! . 

n=l 
-1/2 

71=1 

-1/4 
n=l/4 




exp 



71 
2 



TT 

48 



ln^2 



In^ TT 
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Here, 7 = 0.577215 . . ., 71 = .072815 . . ., and G = .91596 ... are the Euler- 
Mascheroni, Stieltjes, and Catalan constants, respectively. Again, these formu- 
las have classical limit representations looking like equation which we do 
not write down here explicitly. 



4.3 Perspective: A Series by Gosper. 

The paper "On some strange summation formulas" [6] contains some formulas 
like (f25|) below. There might possibly be very short proofs for all these iden- 
tities using fractional sums. The only problem is that there is one single step 
(indicated by the question mark) which we are unable to justify: it is basically 
an interchange of a fractional sum and an infinite series. 

Nevertheless, we give this calculation as a speculation, just to show that it 
is tempting to have a closer look at what else might still be possible. 

Speculation (A Series by Gosper). 
For every b G M, we have the identity 



S(b) = V (-1)" sm^b^ + 7T^{n + l/2y ^ 7rsin6 

"Proof." We start by writing the aforementioned series as a fractional sum 
with —1/2 terms. By plugging in the definition, one easily confirms that 



We will now use the basic identity 

^ = (27) 

v—x 

for every a; e C and n G N, which can be shown in two different ways: The first 
possibility is to sec that for every x S — N, 

^ ^ ^2„+l ^ ^ ^)2„+l + . . . _ i)2„+l ^ (_^)2n+l ^ 

v—x 

since even and odd terms cancel each other. By continued summation and 
Proposition ll.il y'^n+\ _ y'in+\ — ^^"J i/2n+i jg ^ polynomial in 

a:, so equation (1271) must be valid for every a; e C. 

A second way is to use the mirror series from (|20p to calculate 



For polynomials, left and right sum coincide trivially, so (j27p follows immedi- 
ately. 
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Going back to the fractional sum in the odd function 



^^J_sin^^^ (29) 

is holomorphic in the entire complex plane, except for a pole at n = 0, so we 
can develop it into a power series. We get 

-3/4 .^^^ 

S'(6) = - E ( "^""^ + + ^^^"^ + '^5"^ + ' 

n=3/4 ^ 

The next step is critical: we apply the fractional sum term-by-term. Unfortu- 
nately, it is not clear that this manipulation is justified. 

. , -3/4 -3/4 -3/4 \ 

Sib) = -['^ E ""'+^1 E « + ^3 E n^ + ...]. (30) 

ri=3/4 n=3/4 n=3/4 / 

Equations (HH) and (H?]) yield 

sin 6 /3 \ 7rsin6 
5(&) = -^^cot ^-n)^^. 

□ 

This method only works for a certain class of functions which obviously con- 
tains f{n) from ([29]) and other functions like e''"/n'^, but which does not contain 

_ 2 

other simple functions like e . It is an open question to give sufficient condi- 
tions for the validity of this method, i.e., for justification of termwise fractional 
summation as in equation (pO| . 
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